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An adaptive mesh for singular perturbation problems in two 
and three dimensions is reported. The adaptive mesh is generated 
by the solution of potential equations which are derived by 
minimizing the integral I, written, 

I = j[( <V 5 ) 2 +(V n ) 2 }+A v {wJ} + X 0 { (V c *V Ti ) 2 }]dV (1) 

D 

where x(£,rj) ,y (£,ri) represent a mapping from a parameter space 
P, 0<£<I,0<r|<J, where w(x,y) > 0 is given, A y and A Q 
are nonnegative constants, and J, the Jacobian, is written, 


j = liiLiXl 
3(5, n) 


( 2 ) 


In the usual manner, the mesh is constructed by joining points 
in (x,y) corresponding to integer values of £ and n by straight 
lines to form a net of arbitrarily shaped, quadrilateral cells 
( 1 ) • 


The variational formulation is equivalent to Winslow’s 
method (2) when A Q and ^y are zero. The Euler equations are 
those given by Winslow, and their solution maximizes the smooth- 
ness of the mapping. The additional terms modify other attributes 
of the mapping in a similar way. When A Q > 0 , the mapping is 
modified to make it more orthogonal. When Ay > 0 , the mapping 
is modified to make wJ^more nearly constant over the mesh. By 
choosing w(x,y) appropriately, and with Ay , A 0 > 0 , the zone 
size variation and skewness can be controlled. 

In singular perturbation problems, control of zone size 
variation can affect the effort required to obtain accurate, 
numerical solutions of finite difference equations. Consider 
a simple, difference approximation. 


< 



£ i+l - f i-l 
x i+l " x i-l 


(3) 


where i corresponds to £ , 0 _< i < I . In the usual manner, 
the truncation error is written, 
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e = I d£f + 1 dff , )2 

2 dx 2 6 dx 3 ^ "• * 

When f is sufficiently smooth, e is least when x^ = 0 


when f is given, for example, 


( 4 ) 

However, 


f = (1 + exp (x/<5 ) ) 


-1 


(5) 


for which d n f/dx n = 0(6 n ) is not finite for 6=0, the error 
e is bounded only if (x^/6) < 1 in the interval -6 < x < 6. 

An equally spaced mesh with xr sufficiently small satisfies this 
requirement, but one with x^ < 6 only when -6 < x < 6 , and 
larger everywhere else satisfies it with fewer mesh points. 

With w(x) given by, 


w(x) 


1 df 
f dx 


( 6 ) 


minimizing the integral in Eq. 1 causes the mesh spacing to 
approach that given by the equation, 


1 df 

1 


x^ = const. 


(7) 




as Ay increases. As a result, where w is largest the mesh 
spacing is smallest, and vice versa. 

Numerical results for a singular perturbation problem in 
two dimensions are shown in the accompanying figures. A steady 
solution to the equation, 

|| + V • (u<f>) = KV 2 <j> , (8) 


is sought for small values of k with u given. Such solutions 
are obtained when the diffusion and convective transport are in 
balance everywhere. 


u = fv* (9) 

When u is given by, 

u = - | (l+exp((r-r o )/ic)) _ 1 (l+exp(-(r-r o )/K) ) 1 , (10) v. 
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the convective transport term is significantly different from 
zero in an annulus of width k with radius r Q . When w(r,0) is 
given by , 


w(r,6) = u*u 


( 11 ) 


the zones of the computation mesh are made smaller where |u| 
is largest as shown in Fig. 1. 

In Fig. 2, the error in the numerical solution of Eq. 8 
on the adaptive mesh as measured by the maximum norm, 


e 

max 


Max I u 

i» j ~ 


<| v*> 


( 12 ) 


is compared with the error on a fixed rectilinear mesh. The 
accuracy obtained by adapting the mesh can be obtained, in 
most cases, only by a threefold refinement of the regular mesh 
in each coordinate direction. 

The adaptive mesh has been used in calculations of resistive 
magnetohydrodynamic flow in two dimensions with the weight 
function , 


w 


(V x B/B ) 1 2 


The results indicate a significant increase in the maximum, 
representable magnetic Reynolds number. The adaptive mesh 
can be applied easily to other fluid flow problems with the 
appropriate choice of weight function. 

The use of the adaptive mesh in time dependent flow problems 
will be discussed, and results will be presented. 
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Figure 1.- An adaptive mesh with r Q equal to 1/4 and 
k equal to 1/40 the mesh width. The cells are 
concentrated in the region of maximum gradient. 



Average zone size( in inverse mesh widths) 


Figure 2.- Similar scaling of the maximum relative error 
with zone size obtained for both adaptive meshes (like 
the one shown in Fig. 1) and fixed meshes, but many 
fewer zones are required with an adaptive mesh. 



